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In the limit of perfect nesting, the physics of iron-pnictides is governed by the density wave 
formation at the zone-edge vector M. At high energies, various spin- (SDW), charge- (CDW), 
orbital/pocket- (PDW) density waves, and their linear combinations, all appear equally likely, unified 
within the unitary order parameter of U(4) x t/(4) symmetry. Nesting imperfections and low-energy 
interactions reduce this symmetry to that of real materials. Nevertheless, the generic ground state 
preserves a distinct signature of its highly symmetric origins: a SDW along one axis of the iron 
lattice is predicted to coexist with a perpendicular PDW, accompanied by weak charge currents. 
This "hidden" order induces the structural transition in our theory, naturally insures T s > Tjv, and 
leads to orbital ferromagnetism and other observable consequences. 



The discovery of high-temperature superconductivity 
(HTS) in iron-pnictides [H |2] has sparked an intense ac- 
tivity [3]. Like the cuprates, the pnictides are layered 
systems and exhibit anti-ferromagnetism (AF) at zero 
doping (x = 0), followed by HTS beyond some finite x 
[31 3] . Magnetic order in parent compounds consists of 
an AF spin chain along the wave vector (tt, 0) or (0, n) in 
the unfolded Brilliouin zone (UBZ) and an FM spin chain 
along the perpendicular direction 5] . The dynamical ori- 
gin of this AF state is hotly debated: Within the itinerant 
electron model, the magnetic transition is ascribed to the 
SDW instability, enhanced by the near-nesting among 
electron and hole pockets of the Fermi surface (FS) [BJ- 
[9]. To insure "striped" spin order, only one electron 
pocket is involved in SDW, and the spin- wave anisotropy 
arises from the electron pockets' finite ellipticity [TUllllj . 
In contrast, within the localized Heisenberg-type model 
[T2J [T3] various frustrated couplings Ji a , Ju, Ji between 
neighboring spins conspire to produce the observed mag- 
netic order and the magnon anisotropy [HI [15] . 

In addition, the tetrahedral-to-orthorhombic struc- 
tural transformation is observed, accompanied by the AF 
transition [16l [17] . The AF ordered moment is linearly 
proportional to orthorhombicity upon change in x, and 
both transitions disappear for x > x c [18 ]. Magnetoelas- 
tic coupling was suggested as being responsible for the 
close relation between two transitions [E5]. In this ap- 
proach, the structural transition is driven by magnetic 
interactions [20]. However, in the 1111 compounds, the 
structural transition temperature T s is consistently above 
the AF one, Tjy, at any x [5]. Furthermore, the in-plane 
resistivity anisotropy develops well above Tjv in presence 
of uniaxial pressure, and hints at the appearance of a new 
form of order near T S [2T]. One possible explanation for 
T s > Tjv is that magnetic fluctuations are much stronger 
than those associated with structural order. 

In this Letter, we advance another physical picture to 
account for this evident close relation between the struc- 
tural and magnetic transitions: the two are just differ- 
ent faucets of one and the same type of ordering of much 



higher, J7(4)xt/(4) symmetry. This high symmetry char- 
acterizes the dynamics of pnictides within the high-energy 
regime, extending from the energies of order of the effec- 
tive bandwidth D down to those set by T s ~ Tjy. This 
regime is governed by "perfect" nesting and the tendency 
toward formation of a valley-density wave (VDW) at the 
nesting vector Q, with all of its different reincarnations - 
various spin-, charge-, and orbital/pocket-density waves, 
SDW, CDW, PDW, respectively, as well as their mutu- 
ally orthogonal linear combinations - unified within a 
unitary C/(4) x E7(4) order parameter [25]. At yet lower 
energies, however, as the U(4) x t/(4) symmetry-breaking 
interactions and the deviations from perfect nesting come 
into play, the symmetry is reduced down to that of 
real materials. Nevertheless, provided there is a signif- 
icant segregation of scales in the effective Hamiltonian 
of iron-pnictides between the high-energy 17(4) x U(4)- 
symmetric and the low-energy symmetry-breaking terms, 
the ground-state and its excitations bear a distinct sig- 
nature of their highly-symmetric origin. 

Our picture is based on the itinerant model, and re- 
lies on the hierarchy of energy scales that separate the 
"flavor" -conserving from the "flavor" -changing interac- 
tions of quasiparticles on the FS, composed of two hole 
(hi, h%) and two electron (ei, e-i) pockets (or valleys) 
(Fig. 1). This hierarchy is further assisted by the differ- 
ences in area and shape of different pockets being much 
smaller than their common overall features; hence the 
[/(4) e x U(4)h symmetry. Such hierarchy, quantified in 
[22], does not reflect a deep underlying principle; rather, 
it is an accident of the particular semimetallic character 
of pnictides and a screened Coulomb repulsion [23) . But 
be that as it may, the hierarchy is well obeyed in all par- 
ent compounds and we use it as an organizing framework 
to derive the following results: i) The ground state of par- 
ent pnictides is the combination of a SDW along the wave 
vector (tt, 0) or (0, tt) in the UBZ and a spin-singlet den- 
sity wave (DW) along the perpendicular direction; ii) The 
spin-singlet DW is predominantly a PDW, with a tiny ad- 
mixture of a CDW, and is imaginary, i.e. it represents a 
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modulated pattern of weak currents on inter-iron bonds. 
This PDW is difficult to detect and is dubbed the "hid- 
den" order; Hi) The imaginary PDW at Q = (it, 0) (or 
(0, 7r)) induces real CDW at 2Q = (0,0), different from 
the CDW similarly generated by the SDW. The result- 
ing broken orbital symmetry between e x and e y pockets 
(Fig. [IJ drives the observed tetragonal-to-orthorhombic 
transition; and iv) The predicted electronic structure of 
the ground state has numerous observable consequences, 
some of which we explore. Our results are generic for the 
1111 and 122 materials, and - with details changing from 
one compound to another - the overall physical picture 
should be universally applicable. 

M 2 = (0, n) 



V = (0,0) 
ffhl 



Mi = (tt,0) 



FIG. 1. Fermi pockets in the UBZ of iron-pnictides. Two hole 
pockets h\ and /12 are centered at the F = (0,0) point. The 
electron pockets e x and e y are centered at the nesting vectors 
Mi = (?r, 0) and M2 = (0,7r), respectively. The hi, e x , and 
e y pockets are assumed to be perfectly nested to the leading 
order, while /12 is larger than these three; this difference, how- 
ever, is small compared to the overall bandwidth D, as is the 
finite but small ellipticity of e x and e y pockets [7lll0|. 



First, we set up the problem: the band structure can be 
described by the five 3dFe and three pPn orbitals tight- 
binding model [7J, resulting in the FS of Fig. [T] Our 
point of departure is the Hamiltonian H = Hq + H\y : 
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where a, a, and (3 are the spin, hole (h) and (e) pocket 
indices, respectively (Fig. [TJ (3 = x, y for e bands, a = 1, 2 
for h bands) and and are the density operators 
within the e and h pockets [22] . 

H ([I]) describes the high-energy physics of pnictides. It 
contains only the density- density, flavor-conserving inter- 
actions between different pockets, W < D [24]. In con- 
trast, the flavor-changing interactions and the variations 
among Ws in different pockets are all <C D, as long as 
the Hund coupling Jh <C Ud, the Hubbard repulsion on 
d-orbitals [7j|22]. Furthermore, we also initially assume 
perfect nesting, i.e., e£ = e| = -e£ +Ml = -4+m 2 = £ k, 
since the differences among h and e bands are also <C D. 

H (JlJ has a large Z7(4) e x L7(4)/j symmetry, made man- 
ifest by introducing annihilation operators Cn and d v 



to represent h and e pockets, respectively, with fi. v 
1, . . . , 4 labeling both spin and band indices: 
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The interaction H\y drives a VDW formation at the nest- 
ing vectors Mi and M2 (Fig. 1). The order parameter is 
a 4 x 4 matrix A M „ whose 16 complex elements describe 
various SDWs, CDWs, PDWs, and their linear combina- 
tions that gap the FS below some temperature Ty: 
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Integrating out the fermions yields an effective action 
Sa for bosonic fields Ap„. S^A^] has the f/(4) e x 
U(i)h symmetry, spontaneously broken at Ty. Near Ty, 
a Ginzburg-Landau (GL) expansion in A M „ gives |25j: 
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where iV(0) is the density of states of a Fermi pocket and 
{uj n } are Matsubara frequencies. For T < Ty, a < and 
F has a nontrivial minimum for A' A = A 2 , = —a/ ft. 
The solution is A = A U, where U is a 4 x 4 unitary ma- 
trix. At this stage, the four complex 4- vectors comprising 
U describe a plethora of SDWs, CDWs, PDWs, etc., and 
all their mutually orthogonal linear combinations. 

Now, we are ready to confront the real iron-pnictides. 
We turn on all low-energy (<C D ~ W) features ignored 
in ([T]) - differences among Ws, flavor-changing vertices, 
nesting imperfections, and the like [M] - and proceed to 
systematically decode their effect on the £7(4) e x U{A)h 
symmetric theory [26] . The most important among these 
is the interband vertex G 2 , which generates the s su- 
perconductivity as the nesting subsides [SJ [7J : 

Gf 1 c t 1 cJ(d 2 di + d 4 d 3 ) + G e 2 h2 clc\{d 2 d l + d 4 d 3 )l + h.c. 



where G e 2 ha = G e / 1 



27] . The leading order 
correction AF G2 (Fig. [2ja)) to the free energy F ([2J is 

~ n(0) 2 {Gf 1 (AnA 22 + Ai 3 A 24 - Ai 2 A 2 i - Ai 4 A 23 ) 

+ Gf 2 (A 3 iA 42 + A33A44 - A 32 A 4 i - A34A43) + h.cX , 
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FIG. 2. (Color) The leading order corrections to F due to 
(a) e — h and (b) e — e interactions, (c) The ground state 
of parent iron-pnictides. The red and black arrows depict 
iron spins and the [±J7", ±|»7] current pattern, respectively. 
The ground state combines two orders: SDW along (tt, 0) and 
the modulated current DW at the wavevector (0,7r), i.e., the 
"hidden" order. 



11(0) « iV(0)ln(f) HZ]- The Cauchy inequality man- 
dates AF G2 > -n(0) 2 A§ (\Gl hl \ + \G? 2 
ity holds when: a) for Gf 1 > 0, A 22 = 



^12; A 2 4 



-A 13 , A 23 



The equal- 

&u, A 2 i = 
A^ 4 . The VDW involving the 
h\ pocket is then the mixture of real SDW and imaginary 
spin- singlet DW; b) for Gf 1 < 0, the VDW involving hi 
pocket is similarly the mixture of imaginary SDW and 
real spin- singlet DW. The same holds for the h 2 pocket. 

Consequently, G 2 fixes the phases of different DWs. 
One expects that both G% , G e 2 h2 > 0, as the prerequisite 
for high T c s H superconductivity. Hence, the ground 
state of parent compounds must be composed of cither 
real SDW(s) or imaginary spin-singlet DW(s); the lat- 
ter is a general combination of PDW and CDW, in the 
nomenclature of [25]. The real DWs are oc cos(M ■ r), 
with peaks and troughs on the iron sites (Fig. |2jc)). 
In contrast, the imaginary spin-singlet DW breaks time- 
reversal and lattice translation symmetries along M, 
leading to charge/orbital current DW on iron bonds. 

But which one is it, spin-triplet (SDW) or spin- 
singlet (PDW/CDW) density wave? We must con- 
sider next the flavor-changing p-h analogue of G 2 , G\\ 

£a,0=l Gl^L ( e ^ e L' + e V<r e la>) h P<?'- G l (< G 2 

[22] ) generates AF Gl from a diagram similar to Fig.[2ja): 
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Here, it is useful to introduce 2 x 2 Gi matrix 
/ Gl 1 IZeG] 2 ^ 



\TZeG 21 G\ 



22 



Since the phases of DWs are fixed by G 2 , only the real 
parts of G"^ contribute to F. Hence, Gi is real and sym- 
metric, and has two real eigenvalues Ai, A 2 , with associ- 
ated real eigenvectors V\ and v 2 . From Q, AF Gl = for 



SDW and is minimized for the state composed of: a) if 
Ai, A 2 > 0, two real SDWs; b) if Ai, A 2 < 0, two imaginary 
spin-singlet DWs, with AF Gl = (Ai + A 2 )II(0) 2 A§ < 0; 
and, c) if Ai < and A 2 > 0, one real SDW and one imag- 
inary spin-singlet DW, with AF Gl = Ain(0) 2 Ag < 0. 
Experimentally, there is only a single SDW at (it, 0) . 
This implies option c): with majority of GiS rather small 
[22], this is to be expected, once we include the (weak) 
electron-phonon coupling [3J 3] and large polarizability 
of pnictide p orbitals [7J. In this case, the leading order 
contribution of AF° 2 and AF Gl to F is 



F 



a(A 2 DW + A| SDW ) + A 1 n(0) 2 A| SDW . 

Asdw aud iAssdw describe the SDW and the imaginary 
spin-singlet DW, respectively, while 

a(Ts D w) - a(TssDw) + AiH(O) 2 
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set the corresponding transition temperatures. As long 
as |A X | -C W [22], Tssdw ^ ?sdw and A SSDW ^ A SDW . 

Consider now v\ — (a, 6), the (real) eigenvector asso- 
ciated with Ai < 0. AF Gl is minimized by 



A(0) = A 



ial —bo~ r 
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where a n = a ■ n, h is an arbitrary unit vector reflecting 
the SU(2) spin symmetry of our theory, and is an ar- 
bitrary angle, signaling an additional degeneracy in the 
Hamiltonian. The second matrix in |5| is a rotation by 
which mixes e x and e y pockets: 

ei = cos0ez — sin0ej, , e 2 = cos0e y + sin0ea; . (6) 

In the state described by ([5]), e\ and e 2 couple to h n — 
ah\ + bh 2 and h p = ah 2 — bhx, respectively, to form two 
DWs. Finally, this remaining 0-degeneracy is lifted by 
the density-density repulsion between e x and e y pockets: 



(did! 



dld 3 



d\di 



with W£ > 0. The leading order contribution to F, 
AF Wk , follows from Fig. [2^£>) , and contains two fermion 
loops, each with three legs. Were the nesting perfect, the 
loop integral would be independent of leg indices, and, 
upon summation over hole indices, the contribution of 
each loop would be oc A^A, but still independent of 0. 

In real pnictides, however, the outer pocket h 2 deviates 
significantly from hi and perfect nesting (Fig.[T|) 3ilTII25]. 
To account for this, we set e£ 2 = e^ 1 + 77, 77 < W < D. 
At the leading order in 77, the 0-dependent term of each 
fermion loop in Fig. p[6) is now finite and contributes 
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u,k 



1 



v 



= 2p v -> 
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[(asin0) 2 +(6cos0) 2 ] oc (ab) 2 + {b 2 



a 2 ) 2 cos 2 sin 2 . 
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Since generally \a\ ^ \b\, AF Wk is minimized for 6 = 
or it/2. Thus, the preferred ground state combines a real 
SDW in one direction and an imaginary spin-singlet DW 
along the perpendicular direction. The nature of this 
imaginary spin-singlet DW depends on the form of v\ 
[3D]. If a s» —6, the spin-singlet DW is predominantly a 
PDW, translating into a purely orbital current pattern. 
However, unless a — — 6, there is also an accompanying 
charge current DW, depicted in Fig. [2jc) . This current 
DW can be weak for generic a ~ —b but should be observ- 
able and is the main prediction of this Letter. Since the 
charge current DW interacts with the underlying lattice 
more strongly than the pure PDW, it favors an additional 
modulated structural pattern along (0, it), on top of the 
one tied to the SDW along (ir, 0). The apparent absence 
of such pattern in pnictides suggests that indeed a w — b 
and the PDW dominates the imaginary spin-singlet DW. 

With two DWs present at Q = Mi and M2, a real 
CDW at 2Q = (0, 0) is induced as a next harmonic [51] . 
First, this is illustrated within a two-band model, with 
one h and one e pocket. In the mean- field approximation: 



Hmf - 
F < Fmf 



Hmf)mf ■ 



U(fi 2 h +n 2 e ), 



(7) 



U is the intrapocket repulsion and <5£ is the relative shift 
of h and e self-energies. Here we assume the e pocket 
dispersion is £k = k 2 /2m — eg. For SE <C Aq, 



Fmf = a(5£)Tr(A t A) + 0(A 4 ), 



(5n e ) = -(Sn h ) « N(0)5E, 

(Hu ~ Hmf)mf w2JV(0) (<5S) 2 (1 + N(O)U) . 

The r.h.s. of Q is minimized when SE = A§/(8e (l + 
N(0)U)), and hence, (5n e ) = -(Sn h ) = g^q^g^j - 

In a realistic four band model, with the induced CDWs 
at 2Q, a lengthy but straightforward algebra yields [2S] 
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As shown earlier, Asdw < Apdw, and thus (e y e y ) > 
(e x e x ) implying unequal occupancy and splitting of d xz 
and d yz orbitals. Consequently, the induced real CDW 
at 2Q = (0, 0) is an orbital ferromagnet which breaks 
the C4 symmetry while preserving the lattice transla- 
tion symmetry, and can be naturally identified as the 
source of the observed tetragonal-to-orthorhombic distor- 
tion. Since the CDW arises simultaneously with the mod- 
ulated DWs, T s = T PDW > T N = T SDW . For T < T SDW 
and < x <C x c , Eqs. ^ also result in orthorhombicity 



cx AgDwKaO; in agreement with [T5]. Additional support 
for this picture of structural deformation comes from the 
universal scaling of magnetization [32) . 

In summary, we have shown that the high-energy 
c7(4) x U{A) symmetry in iron-pnictides naturally leads 
to the prediction of a "hidden" orbital current DW order 
in parent compounds near the Ax = quantum critical 
point and have explored some of the observable conse- 
quences. 
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